Chapter 9 

Nuclear Reactions 

Michael WiescheiQand Thomas Rauschei0 



Nuclear reaction rates determine the abundances of isotopes in stellar burning pro- 
cesses. A multitude of reactions determine the reaction flow pattern which is de- 
scribed in terms of reaction network simulations. The reaction rates are determined 
by laboratory experiments supplemented by nuclear reaction and structure theory. 
We will discuss the experimental approach as well as the theoretical tools for ob- 
taining the stellar reaction rates. A detailed analysis of a reaction is only possible for 
a few selected cases which will be highlighted in this section. The bulk of nuclear 
reaction processes is however described in terms of a statistical model approach, 
which relies on global nuclear structure and reaction parameters such as level den- 
sity and mass and barrier penetration, respectively. We will discuss a variety of 
experimental facilities and techniques used in the field, this includes low energy sta- 
ble beam experiments, measurements at radioactive beam accelerators, and neutron 
beam facilities. 



9.1 Nuclear Reactions in Astrophysical Environments 

Nuclear reactions are the engine of stellar evolution and determine the overall pro- 
duction of the long-lived radioactive isotopes in a variety of nucleosynthesis pat- 
terns. A detailed understanding of the characteristic production and depletion rates 
within the framework of the nucleosynthesis process is crucial for reliable model 
predictions and the interpretation of the observed abundances. 

There are several experimental and theoretical challenges in obtaining stellar re- 
action rates. The interaction energies in stellar environments extend from basically 
zero projectile energy up to only several MeV. This is especially challenging for the 
measurement of the relevant reaction cross sections which can be extremely small, 
especially for reactions with charged projectiles. This also makes theoretical pre- 
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dictions extremely difficult because several reaction mechanisms (see below) may 
compete and simpler approximations may only be of limited use. 

Another challenge arises from the fact that nuclear burning at high temperature 
produces very short-lived isotopes which subsequently decay to long-lived and sta- 
ble species. Current experimental technology can only access a fraction of those yet 
and is still limited in obtaining detailed information on their properties. The possi- 
bility to measure cross sections of reactions with such highly unstable nuclei is even 
more limited currently. Thus, investigations of nucleosynthesis in high-temperature 
environments largely rely on theoretical models, which not only have to treat the 
reaction mechanisms properly but also are required to predict nuclear properties far 
from stability. 

Moreover, depending on the actual plasma conditions, reactions in an astrophysi- 
cal plasma may proceed fundamentally differently from those in the laboratory. This 
is due to two effects. On one hand, laboratory nuclei are always part of atoms or 
molecules whereas astrophysical nuclear burning involves fully ionized nuclei im- 
mersed in a cloud of free electrons (and photons). The Coulomb charge of a nucleus 
is partially shielded by the surrounding electrons but this shielding (or screening) 
effect will be different for an atom or molecule and a plasma because of the differ- 
ent electron distribution and kinetics. While theoretical cross sections always imply 
bare nuclei, the values have to be appropriately converted (also based on a theo- 
retical treatment of different screening mechanisms) for comparison to low-energy 
laboratory cross sections and for application in astrophysical plasmas. Addition- 
ally, the quantum mechanically and geometrically different electron distribution in 
a plasma directly affects electron capture reactions. For example, nuclei such as 7 Be 
or ^Ti, decaying by capturing an electron from the K-shell of the atom under lab- 
oratory conditions will not be able to do this in a stellar plasma. Instead, electron 
capture inside a star involves capturing a free electron from the plasma, which is 
more unlikely and therefore the terrestrial half-life can be shorter than the one in a 
stellar environment (Iliadis 2007 , |Johnson et aT]|1992[ ). 

Finally, due to the high photon and matter densities in astrophysical environ- 
ments, nuclei very quickly reach thermal equilibrium with the surroundings by exci- 
tation and de-excitation via photons and by collisions. In most cases, this happens on 
a shorter timescale than that of nuclear transformations (one of the exceptions being 
isomeric states). Consequently, the nuclei involved in the reactions occur not only 
in their ground states, as in the laboratory, but also their excited states are populated 
with a probability involving the Boltzmann factor. So far, reactions on excited states 
can only be treated theoretically. The population of excited nuclear states does not 
only depend on the plasma temperature but also on the structure of a nucleus. Nuclei 
with isolated levels reachable within a few keV projectile energy will exhibit pro- 
nounced thermal population already at low plasma temperature. This is especially 
important for modern s-process studies which require high accuracy knowledge of 
neutron capture rates. Thermal effects are also important in decays and neutrino re- 
actions because the available phase space of the reaction products is altered, leading 
to a modification of the rate. For example, electron capture rates in the stellar core 
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collapse are enhanced at temperatures T > 1 .5 Me V because of the unblocking of 
low-lying neutron states by thermal excitation (Cooperstein et al. 1984). 



9.1.1 Reaction Networks and Thermonuclear Reaction Rates 

The change of abundances Y with time due to nuclear processes is traced by cou- 
pled differential equations. To be fully solvable, the number of equations N must 
equal the number of involved nuclei acting as reaction partners and thus an equation 
matrix of size N 2 has to be solved. Such a set of coupled equations is called reaction 
network and can generally be written as 

p .y/<< p ^p.'VE;^''rI;'' ;i r,....J , (9.1) 

where 1 < i < N numbers the nucleus, ;Ay is the jth rate for destruction or cre- 
ation of the ;th nucleus without a nuclear projectile involved (this includes sponta- 
neous decay, lepton capture, photodisintegration), and ,-r, is the rate of the y'th re- 
action involving a nuclear projectile and creating or destroying nucleus i. Similarly, 
we have three-body reactions where nucleus i is produced or destroyed together 
with two other (or similar) nuclei. Reactions with more participants (denoted by . . . 
above) are unlikely to occur at astrophysical conditions and are usually neglected. 
The quantities jPj, fPj, and \Pjk are positive or negative integer numbers specifying 
the amount of nuclei i produced or destroyed, respectively, in the given process. As 
shown below, the rates A, r, and r contain the abundances of the interacting nuclei. 
Rates of type A depend on one abundance (or number density), rates r depend on 
the abundances of two species, and rates r on three. 

Using abundances Y instead of number densities n = YpN\ (where p is the 
plasma density) has the advantage that a change in the number of nuclei in a given 
volume due to density fluctuations is factored out and only changes by nuclear pro- 
cesses are considered. Using abundance changes, the total energy generation rate 
per mass due to nuclear reactions can easily be expressed as 

e = -^F,A A M ! c 2 , (9.2) 

i 

with the rest masses M,c 2 of the participating nuclei. 



The rates ,A ; - appearing in the first term of Eq. 9.1 are reactions per time and vol- 
ume, and only contain the abundance Yj. For example, ,-Ay is simply niLj = YjpN^Lj 
for j3-decays. The factor Lj = (ln2) /■'Tin is the usual decay constant (with the unit 
1/time) and is related to the half-life ^T^m of the decaying nucleus j. It has to be 
noted that some decays depend on the plasma temperature and thus Lj is not always 
constant, even for decays. 
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Two-body rates r include the abundances of two interacting particles or nuclei. In 
general, target and projectile follow specific thermal momentum distributions dn\ 
and dii2 in an astrophysical plasma. With the resulting relative velocities V] — \i, 
the number of reactions per volume and time, is given by 

r\2 = J ff(|vi -v 2 |)|vi -y 2 \dn x dn 2 , (9.3) 

and involves the reaction cross section a as a function of velocity /energy, the relative 
velocity vi — \2 and the thermodynamic distributions of target and projectile dn\ 
and dn2- The evaluation of this integral depends on the type of particles (fermions, 
bosons) and distributions which are involved. 

However, many two-body reactions can be simplified and effectively expressed 
similarly to one-body reactions, only depending on one abundance (or number den- 
sity). If reaction partner 2 is a photon, the relative velocity is always c and the 
quantities in the integral do not depend on dn\. This simplifies the rate expression 
to 

X x =L y {T)n x , (9.4) 

where Ly(T) stems from an integration over a Planck distribution for photons of 
temperature T. This is similar to the decay rates introduced earlier and therefore we 
replaced r by X in our notation and can include this type of reaction in the first term 



of Eq. 9.1 A similar procedure is used for electron captures by protons and nuclei. 
Because the electron is about 2000 times less massive than a nucleon, the velocity of 
the nucleus is negligible in the center-of-mass system in comparison to the electron 
velocity (|v nU cieus - v e i e ctron| ~ |v e i e ctron|)- The electron capture cross section has to 
be integrated over a Fermi distribution of electrons. The electron capture rates are a 
function of T and n e = Y e pNA, the electron number density. In a neutral, completely 
ionized plasma, the electron abundance Y e is equal to the total proton abundance 
Y e = Y,i Z{Yj and thus 

^nucleus, ec = ^ec P^eV^nucleus • (9-5) 

Again, we have effectively a rate per target L (with unit 1/time) similar to the 
treatment of decays earlier and a rate per volume including the number density 
of only one nucleus. We denote the latter by A and use it in the first term of Eq. 



9.1 This treatment can be applied also to the capture of positrons, being in ther- 
mal equilibrium with photons, electrons, and nuclei. Furthermore, at high densities 
(p > 10 12 gcm~ 3 ) the size of the neutrino scattering cross section on nucleons, nu- 
clei, and electrons ensures that enough scattering events occur to lead to a continu- 
ous neutrino energy distribution. Then also the inverse process to electron capture 
(neutrino capture) can occur as well as other processes like, e.g., inelastic scattering, 
leaving a nucleus in an excited state which can emit nucleons and a particles. Such 
reactions can be expressed similarly to photon and electron captures, integrating 
over the corresponding neutrino distribution. 

In the following, we focus on the case of two interacting nuclei or nucleons as 



these reactions will be extensively discussed in Sees. 9.3 and 9.4 This will result in 
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an actual two-body rate r to be used in the second term of Eg. |9~T| Here, we mention 



in passing that Eq. 9.3 can be generalized to 3 and more interacting nuclear species 



by integrating over the appropriate number of distributions, leading to rates f and 



higher order terms in Eq. 9. 1 



Turning our attention back to two-body reactions, we note that the velocity dis- 
tributions can be replaced by energy distributions. Furthermore, it can be shown 



that the two distributions in Eq. 9.3 can be replaced by a single one in the center- 



of-mass system. This time the resulting expression describes a rate r including two 



abundances (or number densities) and showing up in the second term of Eq. 9. 1 
The rate r is defined as an interaction of two reaction partners with an energy dis- 
tribution <p(E) according to the plasma temperature T and a reaction cross section 
o(E) specifying the probability of the reaction in the plasma: 



1 + <5i2 Jo 



a(E)<j){E)dE . (9.6) 



The factor 1 /(l + 6*12) with the Kronecker symbol 8 is introduced to avoid double 
counting. The nuclear cross section is defined as in standard scattering theory by 

number of reactions target - 1 sec - 1 

(7 = . (9.7) 

flux of incoming projectiles 

However, in an astrophysical plasma, nuclei quickly (on the timescale of nuclear 
reactions and scattering) reach thermal equilibrium with all plasma components. 
This allows thermal excitation of nuclei which follows a Boltzmann law and gives 
rise to the stellar reaction rate 

r* =-^| 1 ^\(2J x + l) ro*(E*)4>(E*)e-^dE*\ 

«1«2 1 

where the sum runs over all excited states x of the target (for simplicity, here we 
assume the projectile, i.e. the second reaction partner, does not have excited states) 
with spin J x and excitation energy E x , The quantity G is the partition function of 
the nucleus. The cross section a x includes the reactions commencing from excited 
state x and they are functions of the energy E x relative to this excited state. Cross 
sections <y = <x T=0 measured in terrestrial laboratories do not include such thermal 
effects. At low temperature (e.g., for the s-process) the stellar enhancement factor 
SEF = r* jr will only differ from unity when there are excited states within a few 
keV above the reaction threshold. At the large temperatures reached in explosive 
burning, thermal enhancement can lead to a considerable deviation from the ground- 
state cross section, see also Rauscher (201 1 ). 

Nuclei in an astrophysical plasma obey a Maxwell-Boltzmann (MB) distribution 



Y,{(2J x +l)£a x (E x )<l>(E x )e-*fdE x } , (9.8) 



0(E) = 0mb an d we obtain finally (Rauscher 201 1 k 
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r = Trf«»>' • (9 - 9) 
1 + 612 

<C7V> * =( ^ )1/2( ^ r3/2 G(7 ; )^{ (2 - /v + 1) ■ 

(9.10) 

Here, jj, denotes the reduced mass of the two-particle system and < av >* is the 
stellar reaction rate per particle pair or reactivity. 

As mentioned above the charge of the reaction partners can be screened. For 
most astrophysical conditions this can be included by introducing a screening factor 



/screen, modifying the above rate for bare nuclei (Iliadis 2007 Salpeter et al. 1969 1 



r x =f m r • (9.11) 

The screening factor is derived from the plasma conditions of the specific stel- 
lar environment. At high densities and low temperatures screening factors can en- 
hance reactions by many orders of magnitude and lead to pycnonuclear ignition 
( |Yakovlev et al.[ |2006| ). However, note that the above factorization is not valid for 
vanishing temperatures when nuclei are trapped in a Coulomb lattice. 

Forward and reverse rates are related. Applying the well-known reciprocity the- 
orem for nuclear transitions (Bla tt, et al.| |199l) and further assuming that the re- 
actands in the entrance channel a as well as the reaction products in the exit chan- 
nel b are instantaneously thermalized (the detailed balance principle), the relation 
dHolmes et al.l [T976| |IUadis l |2007) 

< crv >U= |±^g^ ( < av > U , (9-12) 



1 + S ai a 2 Gi H Gb 2 \Ht 

relating the stellar reverse rate to the stellar forward rate. The latter has the reac- 
tion g-value Q. For captures (forward channel a) and photodisintegrations (reverse 



channel b), Eq. 9.12 transforms to 



I G ai G a2 fn a kT\ 3/2 a 



^-I+^TTISfJ ."-<«>^„ ■ (9.13) 

These expressions will not be valid anymore if any of the involved rates was derived 
from a laboratory cross section. They also imply that the detailed balance assump- 
tion is valid. Detailed balance can be violated in nuclei with long-lived isomeric 
states which are not populated or depopulated during regular reaction timescales. 
For these cases, reactions to separate final states have to be calculated and the 



(de)population of these states by photon transitions followed explicitly (Ward et al. 



2002). 



1980). Important examples for such nuclei are 26 Al and 180 Ta (Rauscher et al. 
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9.1.2 Reaction Equilibria 



It is not always necessary to solve a full reaction network (Eq. 9.1 1 including all 
the rates. On one hand, simplifications can often be made by omitting slow reac- 
tions which will not contribute significantly during the timescale of the astrophys- 
ical event. These are, for example, charged-particle reactions on heavy targets in 
hydrostatic stellar burning. On the other hand, high temperature can establish reac- 
tion equilibria. When both forward and reverse reactions become sufficiently fast 
to reach equilibrium with abundances set at equilibrium values. The equilibrium 
abundances of nucl ei can be derived by using the relations |9.12| and |9.13| in the net- 



work equation 9.1 and assuming 7 = 0. Somewhat depending on the density, for 
T > 4 — 5 GK all reactions are in a full nuclear statistical equilibrium (NSE) and 
the abundances are given by 

^AiYj = 1 ^ZjYi = Y e (9.15) 

( i 

with Z,-, Ni, A,-, and B; being the charge, neutron number, mass number, and the 
binding energy of the nucleus i, respectively, the atomic mass unit m u , and the 
abundances of free neutrons Y n , free protons T p , and free electrons Y e . Here, it is 
assumed that reactions via the strong and electromagnetic interactions are in equi- 
librium while the weak interaction is not. Therefore, Y e can still be time-dependent 
and thus also the resulting NSE abundances Y,. 

At T < 4 GK and/or low densities only some reactions may be in equilibrium 
while others are too slow. This gives rise to the so-called quasi-statistical equilib- 
rium (QSE) where only groups of nuclei are equilibrated and those groups are con- 
nected by slower reactions which are not in equilibrium ( |Hix et aTj |1999|l. Abun- 



dance ratios within a QSE group can be determined by application of Eq. |9.14 while 



the connecting, slow reaction determines the amount of matter in each group relative 
to the other groups at a given time. QSE occurs in low temperature, low density Si- 
burning and in O-burning of stars. Often, the slowest rate falling out of equilibrium 
first is that of the strongly density-dependent triple-a reaction. 

A special case of a QSE is the waiting-point approximation, often used in r- 
process calculations ( |Cowan et al.||199Tj|Arnould et al.| |2007). There, the network 
is reduced to neutron capture reactions and their reverse reactions, and j3~-decay 
(with possible release of neutrons). Assuming equilibrated capture and photodisin- 
tegration, QSE within an isotopic chain is obtained and the relative abundances are 
given by 

Y i+ i G i+1 (At + 1 \ 3/2 / 2nh 2 \ 3/2 gnc ap 



^=n n ^^ \ - L ^-\ eir . (9.16) 

Yi a 2Gi \ At J \m u kTj 

The neutron number density is denoted by « n and the neutron capture 2 _va l ue is 
given by the neutron separation energy in nucleus i+ 1: 2ncap = Sn,i+l- The indices 
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i are ordered by increasing neutron number. The j3~ -decays are much slower and 
not in equilibrium. Synthesis of the next element is delayed until the decay of the 
produced isotopes. Typically, only one or two nuclides have significant abundances 
in such an isotopic QSE chain, hence the name waiting-point approximation. 

The advantage of using equilibria is that the rates - and thus the cross sections - 
do not have to be known explicitly. The resulting abundances are completely deter- 
mined by basic nuclear properties and the conditions in the astrophysical environ- 
ment. 



9.2 Relevant Energy Range of Astrophysical Cross Sections 



In the general calculation of the reaction rate according to Eq. 9. 10 the nuclear cross 
section has to be known. Although the integration limits in Eq. 9.10 run from Zero 
to Infinity, significant contributions to the integral only come from a comparatively 
narrow energy range. This is due to the shape of the MB distribution, showing a 
peak around the energy £mb = kT and quickly approaching very small values both 
towards E = and E S> kT. For a slowly varying cross section (as found, e.g., in 
non-resonant neutron-induced reactions), the relevant energy range is simply given 
by the peak of the distribution, Eq = £mb and its width Aq = Amb- For partial waves 
higher than s-waves, the additional centrifugal barrier introduces a stronger energy 
dependence in the cross section and shifts the relevant range to slightly higher en- 
ergy, i.e. E « 0.172T 9 (£+ 1 /2) MeV and A w 0.1947^1+ 1/2 for partial waves 
> (Rau scher et al.|[l997| [Wagoner etal. 1969| l. Charged-particle cross sections 



exhibit a strong energy dependence at energies close to and below the Coulomb bar- 
rier. They decrease by many orders of magnitude towards lower energy. Using the 
astrophysical S-factor 

S(E) = oEe 2 ^ , (9.17) 

with tj being the Sommerfeld parameter describing the barrier penetrability, most 
of the Coulomb suppression is taken out and S(E) is easier to handle because it 
is varying less with energy than a. Inserting definition |9. 17 1 into Eq. |9.10| shows 
that the penetration factor causes a significant shift of the relevant energy range 
towards higher energy. The resulting energy window (the Gamow window, given by 
the Gamow peak appearing when folding the charged particle cross section with the 
MB distribution) can be approximated by Iliadis ( 2007| l; Rauscher et al. ( 1997| l 



E Q = 0.12204 (z\z\\l) 1/3 T 9 2/3 MeV (9.18) 

4 = 4^MI =0.23682 {Z\z\\i) 1/6 Tg /6 MeV . (9.19) 

Here, Tg is the plasma temperature in GK. The idea of a single, relevant energy 
window is only viable for non-resonant cross sections or reactions with broad res- 
onances. Strong, narrow resonances lead to fragmentation of the peak and split it 
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up in several small energy ranges around the resonance energies, with decreasing 
weight towards higher energy. 

It is important to note that Eq. 9. 18 is not always valid. It is based on the assump- 
tion that the energy dependence of the cross section is mainly determined by the 
penetration of the projectile through the Coulomb barrier. However, the dependence 
is dominated by the one of the smallest width in the entrance or exit channel for 
resonant reactions or smallest averaged width in the case of Hauser-Feshbach com- 
pound reactions (see Sec. 9.3.2 1. This smallest width can also be the one of the exit 
channel, leading to a different maximum in the contribution to the reaction rate in- 
tegral than estimated from Eq. |9.18 This is often the case in capture reactions when 
Jprojectiie 3> fy ( jlliadis 2007 Newton et al. 2007 1. Because of the weak energy 
dependence of the y-width, there would not be a Gamow window. Effectively, how- 
ever, the Gamow window is shifted to energies where i^ ro j e ctiie (which is strongly 
energy dependent) becomes smaller than fj,. Since reaction rates at higher temper- 
ature are determined by cross sections at higher energy, the discrepancy between 
Eq. 9.18 and the true maximum of the integrand is more pronounced at high tem- 
perature than at low temperature. Therefore, the relevant energy range for reactions 
in explosive burning should be derived by a proper inspection of the product of the 
(predicted) cross sections and the MB distribution (see Rauscher ( 2010 1 for details). 
For other charged particle captures in astrophysics, often I^rqjectiie <^ ly due to the 
low interaction energy implied by Emb = kT = Ty/ll .6045 MeV, unless for light 
nuclei (with low Coulomb barrier). Regarding neutron captures, although r n 3> ly 
will apply in most cases (unless very close to the reaction threshold), the shape of 
the integrand is mostly determined by the shape of the MB distribution and obvi- 
ously not by any Coulomb penetration. Therefore, the relevant energy window for 
neutrons can still be estimated from the MB distribution as shown above. 
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Having determined the relevant energy range, the cross sections have to be predicted 
by reaction models or determined experimentally. As previously mentioned, often 
measurements for astrophysics prove difficult due to small cross sections or/and 
unstable nuclei involved. However, even if a measurement is feasible, the result- 
ing cross section has to be corrected for effects of electron screening and thermal 
excitation of the target via theoretical models before being used to compute an as- 
trophysical reaction rate. 

Here, we provide a brief overview of approaches to predict low-energy cross 
sections of reactions involving the strong force. Decays and other reactions via the 
weak force are important but cannot be discussed due to limited space. The reader is 
referred to other sources, e.g. ( |M611er et al]|2003t [Fuller et al.||1982||Vogel| [2006) 
and references therein. We also do not cover fission reactions which are important 
in extremely neutron-rich explosive environments where a r-process could occur 
and reach the region of fissionable nuclei. Current predictions of fission barriers, 
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however, carry large uncertainties. For details, see e .g. (|Cowan et al. 1991[ Arnould| 



eTaE] [20071 jPanov et af| [2005] [Goriely et al.j [20091 |Panov et al.| |26l0] i~and 



references therein. We also only discuss reactions between a nucleus and a nucleon 
or an a -particle as the majority of reactions in astrophysics is of that type. 

The interaction of a particle with a nucleus can excite few or many degrees of 
freedom, i.e. transfer energy to few (or none) or to many of the nucleons constituting 
the target nucleus. In nature, all interaction types are, in principle, possible but often 
only one will be dominating at a given interaction energy but with gradual transitions 
from one type to the other within certain energy intervals. For theory, it is simpler to 
consider extreme, idealized cases. Interdependence and interference effects between 
different reaction mechanisms, even if in principle understood, are very difficult to 
predict and especially so for the required large number of reactions with unstable 
nuclei required in astrophysics. In the following we introduce a selection of relevant 
reaction mechanisms considered in literature. The number of degrees of freedom 
which can be excited depends on the number of states or levels present in the system 
formed by projectile and target (Descouvemont et al. 2006) 1 . Therefore, it is helpful 
to distinguish between compound systems with low and high level densities. 



9.3.1 Resonance and Potential Models 

Low level-density systems exhibit no or only few, isolated resonances in the relevant 
energy range. These involve mostly light nuclei which have few, widely spaced 
excited states within several tens of MeV above the ground state and therefore also 
show only few resonances even when the separation energy of the projectile from 
the compound system is large. A similar situation also occurs for heavier nuclei with 
closed shells or heavier nuclei far off stability and close to the driplines where the 
projectile separation energy becomes very low (e.g. in neutron capture on extremely 
neutron-rich nuclei) and in consequence the compound system is formed at very low 
relative energy. 

In principle, isolated resonances can be included by the Breit-Wigner formula 
( |Blatt,etal.|[T99Tl l 

o* = ^ 2 1±1 F JR (9 20) 

4^(2y, + l)(27 proj + l) (jB _ jBres) 2 + i| 1 ' 

where J and E KS refer to the spin and energy of the resonance, co is the de Broglie 
wavelength, and r m is the total resonance width, including the entrance and exit 
widths r* and plus all other open channels. Note that the widths are energy 
dependent. For a narrow resonance, inserting the above in Eq. |9.10| yields 

,ii i 27 + 1 r//i 
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This gives the reactivity in units of cm 3 s _1 mole _1 when the widths and the reso- 
nance energy E KS are inserted in units of MeV (Iliadis 2007 1. (Note that the above 
equations do not involve stellar cross sections. For a true stellar cross section and 
rate, a thermally weighted sum of target states has to be used, according to Eq. 9.8 ) 
However, tails of resonances with the same J may interfere and there may also be 
interference with a direct component (see below). Therefore, additional interference 
terms may have to be added (see, e.g., (Rauscher et al. 1996| l). Furthermore, lo- 
cation of the resonance and the widths have to be predicted from nuclear structure. 
Currently, this is not possible from first principles (except for the lightest nuclei) 
with the accuracy needed in applications. Therefore, this information usually has to 
come from experiments. 

Instead of Breit-Wigner formulas and interference terms, often the R-matrix 
method ( Lane et al.| |1958] l is used. It is applied to parameterize experimentally 
known cross sections with as few parameters as possible, implicitly accounting for 
resonances and their interference. The R-matrix approach can be used to extrapolate 
the nuclear cross section from existing data to the Gamow range as long as nuclear 
structure information about resonance levels and non-resonant reaction contribu- 
tions are included. 

In addition to possible resonance contributions a direct capture process can oc- 
cur. These are fast, one-step processes in which a captured particle directly enters 
the final state. Typical reaction timescales of direct processes are of the order of 
10~ 22 s whereas compound reactions, distributing the energy among a large number 
of nucleons, take of the order of 10~ 16 s. Direct reactions include transfer processes 
where a particle exchange takes place between projectile and target nucleus, and 
capture processes in which the projectile is being fully captured by the target nu- 
cleus. These two reaction types can be treated in ab initio models, determining the 
cross sections from wave functions obtained by solving the Schrodinger equation 
using effective potentials. 

For transfer reactions often the Distorted Wave Bom Approximation (DWBA) 
(Satchler 1983 Glendenning 2004 1 is used, utilizing optical potentials to compute 
the cross sections from the overlap integral of distorted scattering wave functions 
and the bound state wave function. The DWBA implicitly assumes that elastic scat- 
tering is dominant while non-elastic contributions can be treated perturbatively. 

On the other hand, capture reactions can be calculated with a simple potential 
model, which is a first-order approach involving an electromagnetic operator de- 
scribing the emission of photons due to the dynamics in the movement of electric 
charges ( |Descouvemont et al. 2006} . In the potential model the differential cross 
section is proportional to the matrix element defined by the overlap of the final state 
0^3 of the final nucleus and the initial state composed of the target wave function 
a and a (distorted) scattering wave of the projectile % a . This can be decomposed 
into an overlap function S of the target and the final nucleus and a radial integral 
containing the scattered wave Xa( r )> tne bound state wave function of the projectile 



in the target a +A, and the radial form of the electromagnetic operator ^em (Kim 

[etaLirT987l > 
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~d£2 



| < 0/3 1 ^EM \Za K ) | 2 - ^ 1 a+A (r) ^EM (r)j£ W^f 



(9.22) 



The wave functions a +A.( r ) an d X« (0 are obtained by solving the radial Schrodinger 
equation with appropriate effective potentials. 

Both approaches, DWBA and potential model, require a renormalization of the 
resulting cross section through spectroscopic factors S, describing nuclear structure 
effects by the overlap of initial and final state of the system. These spectroscopic 
factors have to be obtained from nuclear structure models or by comparison with 



experiment ( Satchler 1983 Glendenning, 2004). 

Microscopic reaction models are first principle methods, starting from effective 
nucleon-nucleon interactions and treating all nucleons in a Hamiltonian with exact 
antisymmetrization of the wave functions. Because of this, no artificial distinction 
between direct and resonant contributions has to be made. Unfortunately, such reac- 
tion models are limited to systems of few nucleons. Although the Quantum Monte 
Carlo method (Pieper et al. |200 1 1 > is promising, it is currently limited to A < 10 
and not applicable to continuum states. Cluster models have been often used for 
light systems so far (see Descouvemont et al. 2006 Descouvemont| |2003| and 
references therein). They assume that the nucleons are grouped in clusters and use 
cluster wave functions defined in the shell model and computed with an adapted ef- 
fective nucleon-nucleon force. The Resonating Group model (RGM) and the Gener- 
ator Coordinate Method (CGM) are two equivalent implementations differing in the 
definition of the relative wave function of the clusters ( |Descouvemont et al.| [2006 ; 
[Descouvemont 2003). 



9.3.2 Statistical Model 

In systems with high level density p(J,n,E), individual resonances cannot be re- 
solved anymore and an average over the overlapping resonances can be used instead 
(Fig. |9.1| l. Further assuming that the relative phases are randomly distributed, inter- 
ferences will cancel and a simple sum of Breit-Wigner contributions can be replaced 
by a level-density weighted sum of averaged widths (F) over all spins J and parities 
n ( [Descouvemont et al.| [20061 IGadioli "etaLl[T992l ) 

a\E) ' 
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Fig. 9.1 Applicability of the Hauser-Feshbach model to calculate astrophysical reaction rates for 
neutron-induced reactions: Shown are the stellar temperatures above which the rate can be calcu- 
lated from Hauser-Feshbach cross sections (reprinted from Rauscher etal. 1997 with permission). 



xW(J,n,E c )} , (9.23) 
E c = E +£sep,pro E x , (9.24) 
^fin = E c ^sep.fin ^-.t,fm • (9-25) 

This is called the Hauser-Feshbach or statistical model of compound reactions 
( |Hauser et al.[ |1952| l. Width fluctuation corrections W account for non-statistical 



correlations but are only important close to channel openings (Ericson 1960). The 
separation energy £ se p,pro of the projectile in the compound system determines at 
which energy E c the compound system is formed. The averaged width of the exit 
channel (If,) usually includes a sum over energetically possible final states at energy 
Ex fin or an integral over a level density of the final system when individual states are 
not known or numerous. For capture, compound and final system are identical. The 
averaged widths are related to transmission coefficients 3T = 2np (F). The latter 
are calculated from the solution of a (radial) Schrodinger equation using an optical 
potential. (It is to be noted that these potentials differ from the ones employed for 



low-density systems described in Sec. 9.3.1 ) 



The challenge for nuclear astrophysics lies in the determination of globally appli- 
cable descriptions of low-energy optical potentials as well as level densities, masses 
(determining the separation energies), and spectroscopy (energies, spins, parities) 
of low-lying excited states, to be applied for a large number of nuclei at and far 
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from stability. For details on the different properties and the remaining open prob 
lems in their treatme nt (see, e.g., | Descouvemont et al. 2006| Rauscher etal. 1997 



Arnould et al. 2007 1 Goriely et al. 2008 ; Rauscher 201 1 and references therein). 



For a general discussion of the applicability of the statistical model, see |Rauscher 
|eTaIl ( flW7l|200^ ; |Rauschlr1 ( |26TT 



9.4 Experimental Facilities and Techniques 

The experimental determination or verification of nuclear reaction rates requires a 
large variety of facilities and techniques. This is in particular true if one wants to es- 
tablish experimentally reaction rates associated with the production of long-lived ra- 
dioactive isotopes associated with galactic gamma sources. Nuclear astrophysics re- 
lated experiments include low energy high intensity accelerator measurements with 
stable beams to study charged particle reactions of relevance for quiescent stellar 
burning which may possibly lead to the production of 7 Be, 22 Na and 26 Al. High flux 
neutron beam studies to explore neutron induced reactions for the weak and main 
s-process which can be associated with the production of long-lived radioactive iso- 
topes such as 60 Fe and 98 Tc. Real and virtual photon beams are increasingly used 
for probing nuclear reactions associated with explosive nucleosynthesis events such 
as the p-process but can also be used to probe indirectly neutron capture reactions 
associated with the s-process. Intense radioactive beams are the primary tools for 
exploring nuclear reactions and decay mechanisms far of stability which are ex- 
pected to occur in explosive stellar environments and may lead to the production of 
long-lived radioactive elements such as 18 F, 26 Al, 44 Ti and 56 Ni. 



9.4.1 Low-energy Facilities, Underground Techniques 

Low energy charged particle measurements belong to the most challenging exper- 
iments in nuclear astrophysics. The cross sections need to be measured at the ex- 
tremely low energies associated with the Gamow range of quiescent stellar burning. 
This requires to determine the cross sections of proton capture reactions for hydro- 
gen burning in main sequence stars at energies well below 100 keV. Measurements 
for helium burning in red giant stars need to be explored in the 200 keV to 500 keV 
range and heavy ion fusion reactions in subsequent stellar evolution phases need 
to be measured near 1 MeV to 2 MeV center of mass energy. The cross sections 
are extremely low, typically in the femto-barn range, which requires a long time, 
in excess of days, to accumulate a statistical relevant amount of reaction yield data. 
Typical experimental techniques are summarized in the text book literature (Iliadis , 
2007) and will not be discussed here. 

The critical issue with low cross sections is that the yield of reaction events is low 
in the detectors measuring the characteristic gamma or particle radiation produced. 
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This requires using high efficiency detector material with high resolution to sepa- 
rate the characteristic events from random background events. High beam intensity 
is desired to increase the event rate, however it may also increase beam induced 
background on target impurities and is limited by target stability. 

The second critical issue is the background rate in the detector. There are typ- 
ically three different kind of background, cosmic ray induced background in the 
detector environment, natural long-lived radioactivity in the detector material and 
the surrounding environment, and beam induced background on low Z target impu- 
rities and beam defining slits or apertures. This background must be reduced as far 
as possible to identify reaction events in the spectrum. 

Cosmic ray induced background affects the spectra up to very high energies and 
makes it difficult to extract weak signals. That background is the most important 
to remove. Natural environmental background will be strong in an underground en- 
vironment except for salt mine locations. But the characteristic y lines are mainly 
below 3 MeV and can be shielded locally. Neutron background is more difficult to 
absorb and needs special shielding arrangements. Beam induced background de- 
pends critically on the target as well as the choice and preparation of the target 
material. It is difficult to suppress and may require active shielding procedures. 

This can be done by identifying the event electronically by its particular charac- 
teristics such as coincidence conditions in a particular decay sequence, pulse shape 
or timing conditions and reject the background events which do not fulfill these re- 
quirements. This can lead to active background suppression by up to three orders of 
magnitude ( |Runkle et al. 2005 Couture et al. 2008| >. While this clearly helps in 
many cases a more efficient background reduction is desired. 

The high energy cosmic ray induced background can be most successfully sup- 
pressed by operating the experiments in a deep underground environment where the 
cosmic ray flux is heavily reduced. This was demonstrated with the installation of 
the LUNA accelerator facility at the Gran Sasso deep underground laboratory in 
Italy. The cosmic ray induced background was successfully removed and several 
critical reactions of the pp-chains and the CNO cycles were successfully measures 
in the or near the Gamow energy range (Costant ini et al.||2009| l. As a consequence 
of this successful operation new underground accelerator facilities are being pro- 
posed or planned which would allow to cover reactions over a wider energy range 
than available at LUNA. This is of particular importance for an improved R-matrix 
analysis and extrapolation. Higher energies are also of great relevance for the under- 
ground measurements of a capture reactions and stellar neutron sources in helium 
burning. In particular it will also improve the chances for pursuing heavy ion fusion 
reaction studies towards lower energies. There are presently three major initiatives 
for the construction of new underground accelerator facilities. The proposal to es- 
tablish an underground accelerator facility ELENA at the Boulby salt mine in the 
UK seeks to take advantage of the reduced level of neutron and natural activity in 
a salt environment. The disadvantage will be the reduced depth level compared to 
the Gran Sasso location. The second proposal is for the development of a two ac- 
celerator facility DIANA at the DUSEL underground laboratory at Homestake mine 
in South Dakota. The third proposal in debate is the construction of an accelerator 
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Fig. 9.2 The Dragon facility in Vancouver, Canada, is an example of nuclear experiment facili- 
ties, now aimed at experiments for astrophysically-relevant reactions: Radioactive isotopes can be 
selected and accelerated to form projectiles for such reactions of interest. 



facility in an abandoned train tunnel in the Pyrenees mountains at Canfranc, Spain. 
With these facilities the community hopes to address the new and critical questions 
about stellar reaction cross sections and provide the final answer on the nuclear 
engine of stellar evolution. 

However, it has been demonstrated that alternative inverse kinematics methods 
are a very powerful tool in reducing the background. They are based on the tech- 
nique of using a high intensity heavy ion beam on a hydrogen or helium gas target 
and separate the heavy ion recoil reaction products from the primary beam through a 
high resolution electromagnetic mass separator system from the primary beam. This 
method has been demonstrated to be successful at a number of different separator 
facilities such as DRAGON at TRIUMF, Vancouver ( |Vockenhuber et al.| |20"07] > 
(Fig. 9.2) and ERNA at the Ruhr University Bochum pi Leva efaL]|2009) . 

The detection of the recoiling charged particle has a clear efficiency advantage 
compare to the gamma detection. The possible detection of the gamma rays in coin- 
cidence with the reaction products reduces dramatically the backgrounds. However, 
there are several experimental challenges associated with using recoil separators. At 
the low stellar energies, the energy spread and the angular aperture are much larger 
than the acceptance of any of the cited existing recoil separators. In order to measure 
an absolute cross section the transmission of the recoils should ideally be 100% or 
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exactly known. It is also necessary to know precisely the charge state distribution of 
the recoil products. In addition, since the primary beam intensity is typically many 
orders of magnitude larger than the recoiling reaction products, a large spatial sep- 
aration between the reaction products and the beam is required, which is difficult to 
realize for beams with a large energy spread. Therefore, solar fusion reactions are 
particularly challenging to measure with recoil separators and are typically used for 
higher energies and for the helium or heavy ion burning reactions. 

Dedicated next generation separators for low energy nuclear astrophysics studies 
with stable ion beams coming on line are the ST.GEORGE facility in Notre Dame 
( |Couder et al.| |2008| l and the modified and upgraded ERNA facility at CIRCE 
in Caserta, Italy. Both separators feature large angular and energy acceptance and 
will be equipped with high density gas jet targets, which will ensure a well defined 
interaction region. 



9.4.2 Laboratory Neutron Sources 

Many of the observed or anticipated long lived radioactive isotopes in our galaxy are 
produced by neutron induced nucleosynthesis in the weak or main s-process taking 
place in helium and carbon burning stellar environments. This includes 41 Ca, 60 Fe, 
63 Ni, but also more massive isotopes such as 98 Tc and 99 Tc and possibly numerous 
long lived isomers. 

The study of neutron induced stellar reactions leading to the production of such 
isotopes requires high intensity neutron sources with a well defined energy distribu- 
tion to determine the reaction cross sections at stellar energies of a few keV. Neu- 
trons in that energy range can be produced in several ways. Nuclear reactions such 
as 7 Li(p,n) or 3 H(p,n) with high intensity proton beams provided by low-energy 
particle accelerators offer the possibility of tailoring the neutron spectrum to the 
energy range of interest; this has the advantage of low backgrounds. A particularly 
successful approach is to simulate a quasi-stellar neutron spectra in the laboratory. 
In bombarding thick metallic lithium targets with protons of 1912 keV, the result- 
ing neutrons exhibit a continuous energy distribution with a high-energy cutoff at 
E„ = 106 keV and a maximum emission angle of 60 degrees. The angle-integrated 
spectrum corresponds closely to a Maxwell-Boltzmann distribution for kT = 25 keV 



( Ratynski et al. 1988 1. Hence, the reaction rate measured in such a spectrum yields 
immediately the proper stellar cross section. 

Higher intensities can be achieved via photon production by bombarding heavy- 
metal targets with typically 50-MeV electron beams from linear accelerators. When 
these energetic neutrons are slowed down in a moderator, the resulting spectrum 
contains all energies from thermal energy to nearly the initial energy of the electron 
beam. Since the astrophysical relevant range corresponds to only a narrow window 
in this spectrum, background conditions are more complicated and measurements 
need to be carried out at larger neutron flight paths. In turn, the longer flight paths 
are advantageous for neutron-resonance spectroscopy with high resolution. 
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The most intense keV neutron flux is produced by spallation reactions. The LAN- 
SCE facility at Los Alamos is particularly suited for neutron TOF work due to the 
favorable repetition rate of only 12 Hz (Lisowski et al. 1990), and because the ac- 
cumulation of a number of beam pulses in an external storage ring yields extremely 
intense neutron bursts. Accordingly, excellent signal-to-background ratios can be 
achieved. The n_ToF facility at CERN provides high intensity neutron beam pulses 
with a lower repetition rate of 0.4 Hz (Borcea et al. 2003). This has proved highly 
advantageous for a large number of experimental neutron capture studies along the 
s-process path. 

The experimental methods for measuring (n,y) cross sections fall into two 
groups: TOF techniques based on the detection of the prompt capture y rays and 
activation methods. 

The TOF techniques can be applied in measurements of most stable nuclei but 
require a pulsed neutron source to determine the neutron energy via the flight time 
between target and detector. Capture events in the samples are identified by the 
prompt 7 ray cascade in the product nucleus. 

The best signature for the identification of neutron capture events is the total en- 
ergy of the capture gamma cascade, which corresponds to the binding energy of 
the captured neutron. Hence, accurate measurements of (n,y) cross sections require 
a detector that operates as a calorimeter with good energy resolution and is insen- 
sitive to neutron exposure. In the gamma spectrum of such a detector, all capture 
events would fall in a line at the neutron binding energy (typically between 5 and 
10 MeV), well separated from the gamma-ray backgrounds that are inevitable in 
neutron experiments. Such detectors have been successfully developed at the vari- 
ous laboratories using arrays of An BaF2 scintillator detectors with a large number 
of independent detector modules (lleil et aL] 2001 1 and have emerged as standard 
technology for these kinds of measurements. 

A completely different approach to determining stellar (n,y) rates is activation 
in a quasi-stellar neutron spectrum. Compared with the detection of prompt capture 
gamma rays, this method offers superior sensitivity, which means that much smaller 
samples can be investigated. Since it is also selective with respect to various reaction 
products, samples of natural composition can be studied instead of the expensive en- 
riched samples required by the TOF techniques. However, the activation technique 
is restricted to cases where neutron capture produces an unstable nucleus, and it 
yields the stellar rate only for two thermal energies at kT = 25 and 52 keV. This 
method is however particulary powerful in obtaining cross sections for reactions 
producing long-lived radioactive materials which can be identified by their particu- 
lar decay characteristics and signature. This activation technique has been used for a 
variety of measurements. The technique can be applied to short-lived products with 
half-lives in the millisecond range and allows for cross section measurements with 
uncertainties of a few percent. 



9.4 Experimental Facilities and Techniques 



19 



9.4.3 Accelerator Mass Spectroscopy 



Classical activation techniques require a characteristic decay signal associated with 
the decay pattern or the half-life of the produced radioactive isotope. This can be 
difficult in cases where no characteristic gamma or particle decay pattern exists 
or where the decay analysis of the /3 decay signal is prohibited by high back- 
ground activity. In these cases activation analysis through accelerator mass spec- 
trometry (AMS) offers a powerful tool to measure cross sections through ultra- 
low isotope-ratio determination. The AMS method was successfu lly introduced 



for the study of the neutron-capture cross section of 62 Ni(n,y) 63 Ni (jNassar et al. 



2005), and exte nded to other neutro n and charge-partic le-induced reactions, such as 



2S Mg(p,y) 26 Al ( |Arazi et al.||2006| > and 40 Ca(a,y) 44 Ti ( jNassar etal.||2006 i. 

In these cases samples were either irradiated in a neutron spectrum resembling a 
stellar Maxwell Boltzmann distribution or by charged particles of well known en- 
ergies. After the irradiation the samples must be chemically treated to extract the 
radioactive reaction products. This requires some time and limits AMS activation 
studies to more longer lived isotopes. Since isotopic and isobaric interferences may 
represent a major challenge in AMS measurements of irradiated samples, extensive 
background studies for these isotopes are always necessary prior to the irradiations 
in order to demonstrate that the required sensitivity can be reached. In AMS, neg- 
ative ions are extracted from an ion source which have to pass a low energy mass 
spectrometer prior to entering a tandem accelerator. When passing the stripper, pos- 
itive ions are produced while within this stripping process molecular isobars are de- 
stroyed. One positive charge state is selected with a second (high-energy) mass spec- 
trometer system which is optimized for mass, charge and isobar separation trough 
possible combination of dipole magnet separators, Wien-filters, and more recently 
magnetic gas filled separators for improved isobar separation. With such a system 
the concentration ratio of the radioisotope is determined relative to a stable isotope 
by measuring the number of radionuclides relative to the current of the isotopic ions 
in front of the detector, after adjusting the injector magnet, terminal and Wien-filter 
voltage appropriately. By measuring relative to a standard sample of known isotopic 
ratio, factors like stripping yields and transmissions mostly cancel. 

The difficulties with AMS experiments is in the chemical preparation of the sam- 
ple and the sufficient separation of the extracted radioactive ions from background 
events. While AMS is a widely established method with many applications, the 
analysis of the very limited number of radioactive products from low cross section 
reactions remains challenging. Systematic studies are necessary to reduce possible 
uncertainties. 

Dedicated AMS facilities with an established nuclear astrophysics program are 
the Vienna Environmental Research Accelerator (VERA) (Kut schera et al.||l997[ ), 
the Center for Isotopic Research on Cultural and Environmental Heritage (CIRCE) 
in Caserta/It aly (|Terrasi et ak] |2007| l or the Munich Tandem accelerator facility 
( |Knie etaE] [2000) which is optimized for the analysis of more massive radioactive 
isotopes. A new AMS program is presently being developed utilizing the Notre 
Dame tandem accelerator (Roberts on et al.| [2007 ). 
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9.4.4 Radioactive Beam Techniques 

The development of radioactive accelerated beams for low energy nuclear astro- 
physics experiments has been one of the large challenges of the field. The experi- 
mental study of nuclear reactions and decay processes far of stability is necessary 
for the understanding of explosive nucleosynthesis processes such as the rp-process 
in cataclysmic binary systems or the r-process and p-process in the supernova shock 
front. These processes can in particular contribute to the production of long-lived 
galactic radioactivity by primary reaction or also by secondary decay processes from 
the reaction path towards the line of stability. 

For the purpose of studying the origin of long-lived radioactive isotopes in as- 
trophysical environments radioactive beams are utilized in two ways, for producing 
long-lived targets by implantation for subsequent irradiation with neutron, charged 
particle or possibly intense photon beams or for direct reaction measurement in 
inverse kinematics on light ion target materials. The later approach requires well 
defined mono-energetic and intense radioactive beams and a detection system for 
light or heavy recoil reaction products. 

The main challenge in this approach is to produce a sufficiently high intensity 
of radioactive beams which have to be produced on-line as a secondary reaction 
product. This requires high cross sections for the production process and high pri- 
mary beam intensities. A variety of different approaches has been chosen in the past 
to optimize the production efficiency and maximize the intensity of the radioactive 
beams. A technique developed for small scale facilities is the selection of specific 
nuclear reactions tailored for the on-line production of radioactive beams at opti- 
mum conditions. The secondary particles can be used for subsequent nuclear reac- 
tion studies after blocking and separation from the primary beams ( |Kolata et al.| 
1989). The efficient separation of a suitably high intensity beam of radioactive 
species is the most challenging problem for this approach. 

An alternative approach is the ISOL (Isotope Separation On-Line) technique 
where high energy protons are used to bombard heavy ion targets for producing 
a large number of radioactive species through spallation processes. These isotopes 
diffuse out of the target into an ion source for being charged and re-accelerated for 
secondary beam decay or reaction experiments. The method has been proven to be 
very powerful over the years but is limited to isotopes with lifetimes appreciably 
longer that the time necessary for the diffusion transport and ionization process. 
This can be different for different elements because of the associated chemical pro- 
cesses between the isotopes and the surrounding environment. 

The third approach is based on the use of energetic heavy projectiles bombarding 
light or target nuclei fragmenting on impact. This fragmentation process generates a 
cocktail beam of many radioactive species which move forward with high velocity 
because the initial momentum of the primary particles is maintained. For experi- 
ments with a specific secondary particle, it must be selected by fragment separator 
systems which separate and focus the isotopes by magnetic fields and energy loss 
characteristics in heavy wedge materials. For nuclear astrophysics related experi- 
ments the fast beam particles need to be slowed down by energy loss in gas or solid 
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material and re-accelerated to energies corresponding to the temperatures in the ex- 
plosive stellar scenarios. 

There is a number of laboratories which have focused on nuclear reaction stud- 
ies with radioactive beams. The first fully operating radioactive beam laboratory 
based on the ISOL principle was the coupled cyclotron facility at Louvain la Neuve 
which did a number of successful radioacti ve beam studies of rele vance for investi- 
gating the the production of F in novae ( |De Serville et al.||2009| ). These measure- 
ments were complement ed by measuremen ts at the HRIBF facility at Oak Ridge 
using intense 18 F beams (Chae et al. 2006 1. Both facilities produce the radioactive 
species by nuclear reactions on thin production targets, with the reaction products 
being transported into an ion source for producing and subsequently accelerating 
the secondary beam. The intensity is largely limited by target technology and beam 
transport and re-ionization efficiency. 

The premier ISOL radioactive beam facility is ISAC at TRIUMF Canada. The 
primary 600 MeV proton beam is provided by the TRIUMF cyclotron. The reaction 
products are post-accelerated in an RFQ SC LINAC accelerator combination to en- 
ergies of 0.3 to 3 MeV/u. ISAC has successfully performed a number of radioactive 
beam experiments of relevance for explosive hydrogen and helium burning. Most 
notable a direct study of 21 Na(p,y) 22 Mg in inverse kinematics to probe the produc - 
tion mechanism of 22 Na in Ne nova explosion environments ( |D'Auria et al.||2004[ ). 
The facility also runs a successful program with stable beams which was utilized to 
investigate the production of 44 Ti ( Vockenhuber et al. 2007 1. Presently a number 
of studies associated with the production of the long-lived y emitter 26 Al are being 
performed. 

Other ISOL based radioactive beam facilities such as Spiral facility at GANIL 
in Caen, France or REX-ISOLDE at CERN have been used to perform interesting 
experiments for nuclear astrophysics but have been less concerned with the question 
of nuclear production mechanisms for long lived cosmic gamma emitters. 

There have been a number of fast radioactive beam facilities with scientific pro- 
grams in nuclear astrophysics primarily aimed at the study of nucleosynthesis pro- 
cesses far off stability. However the rapid new developments in fast beam physics 
promises a number of new experimental opportunities which can provide benefits 
for studying reactions associated with the production of long-lived gamma emitters 
in explosive nucleosynthesis events. 

There are currently four major fragmentation facilities in the world: GANIL and 
GSI in Europe, NSCL/MSU in the US and RIKEN in Japan. They are all based on 
Heavy Ion accelerators which operate in complementary energy domains. Because 
of the high energy of the fragment products low energy reaction experiments for 
nuclear astrophysics are not possible but the development of indirect techniques to 
determine critical reaction or decay parameters has been the primary goal. In the 
context of long-lived isotopes of astrophysical interest a major contribution was the 
development of fast beams such as 8 B at NSCL/MSU, RIKEN, and GSI for utilizing 
Coulomb dissociation techniques for probing critical reactions such as 7 Be(p,y) 8 B. 
The NSCL and RIKEN also successfully developed a 44 Ti beam for new measure- 
ments of its half-life ( Gorres etal.||1998] l. More half-life measurements of long-lived 
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isotopes such as 60 Fe are presently underway to re-evaluate these critical parame- 
ters. 



9.5 Specific Experiments 

The complexities of the experiments and the uncertainties in the experimental re- 
sults affect the reliability of model predictions on the nucleosynthesis of long-lived 
radioactive species. In particular recent studies of critical nuclear reactions and de- 
cay processes exhibit considerable differences to earlier studies which so far have 
been the reference point for nucleosynthesis simulations and predictions for long- 
lived radioactive isotope abundances in stellar burning environments. It is therefore 
important to carefully evaluate the experimental results and clarify possible discrep- 
ancies and inconsistencies in the data. This section will discuss the present status of 
the experimental reaction rates and evaluate future opportunities to improve the ex- 
isting data base. 



9.5.1 Experiments with Stable Beams 

Many of the long-lived radioactive gamma emitters in our universe have been pro- 
duced by radiative capture reactions on stable isotopes. The best known examples 
are 26 Al, which is primarily formed by proton capture on stable 25 Mg isotopes, 
25 Mg(p,y) 26 Al, and 44 Ti which is most likely produced via alpha capture on sta- 
ble 40 Ca isotopes, 40 Ca(a,y) 44 Ti. Extensive measurements using in-beam y spec- 
troscopy techniques have been made for both reactions and have formed the basis 
for earlier reaction rate compilations. 

The low energy reaction cross section of 25 Mg(p,y) 26 Al is characterized by sev- 
eral resonances with energies between 30 keV and 400 keV . The reaction rate is 
directly correlated to the strengths coy of the resonances. The strengths for the res- 
onances above 190 keV have been determined from the on-resonance thick target 
yield in radiative capture measurements ( |Elix et aL 1979[ Iliadis et al. 1990). The 



strengths of lower energy resonances are estimated on the basis of single particle 
transfer reaction studies. Of particular importance are three resonances at 90 keV, 
130 keV and at 190 keV which determine the reaction rate at temperatures typical 
for stellar hydrogen burning in AGB stars and nova explosions. Because the low en- 
ergy radiative capture measurements have been handicapped by cosmic ray induced 
background, an alternative measurement was done using the AMS technique to an- 
alyze the number of 26 Al reaction products after irradiation at resonance energies 



(Arazi et al. 2006). The experiment was successful and confirmed the resonance 
strengths of the known resonances at 304 keV, 347 keV, and 418 keV resonance 
energy. However the results indicated a substantially lower strength for the critical 
resonance at 190 keV. This would reduce the reaction rate by about a factor of five at 
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the temperature range between 0.2 and 1 .0 GK. This result introduced a large uncer- 
tainty in the reaction rate which motivated a new experimental study at LUNA in the 
Gran Sasso laboratory using in-beam gamma spectroscopy techniques with a variety 
of high efficiency and high resolution gamma detector devices. The measurements 
confirmed earlier gamma spectroscopy studies of the strengths of higher energy res- 
onances (|Elix et al.||1979||riiadis et al.[|1990"| ) tabulated in the NACRE compilation 
(Angulo et al. 1999 1. The new results are being prepared for publication. Parallel 
to the gamma spectroscopy measurement, the irradiated samples were analyzed for 
their 26 Al content using AMS techniques. The AMS measurements were performed 
at the CIRCE facilities. Excellent agreement is demonstrated for the resonance at 
304 keV, additional experiments are being pursued for lower energy resonances to 
address the inconsistencies in the strength determination for the 190 keV resonance. 

The 40 Ca(a, y) 44 Ti reaction is considered to be one of the major production re- 
actions for 44 Ti in supernova shock front nucleosynthesis. The cross section for this 
radiative capture process has been explored in a number of in-beam gamma spec- 



troscopy studies down to center of mass energies of 2.5 MeV ( Simpson et al. 1971 



|Cooperman et al.[ |1977| l. The cross section is characterized by a large number of 
resonances and the initial reaction rate determinations were based on an analysis 
of resonance strengths. Despite the high level density in 44 Ti, it was noted that the 
experimental reaction rate is substantially smaller than the reaction rate based on 
statistical model Hauser Feshbach predictions ( |Rauscher et al. 2000b ). The reac- 
tion was studied independently using a thick He-g as cell target and cou nting the long 



lived 44 Ti reaction products by AMS techniques (jNassar et al.| |2005[) to determine 



the integral yield over an energy range of 1 .7 to 4.2 MeV. The extracted reaction 



rate is substantially higher than the ones discussed in the literature (Rauscher et al. 



2000b ). A more recent study of the reaction using inverse kinematics techniques was 
performed at the ISAC facility at TRIUMF, Vancouver, separating the 44 Ti reaction 
products on-line with the DRAGON recoil separator. The measurements covered 
the energy range of 2.3 MeV to 4.2 MeV (center of mass) in more than 100 small 
energy steps. The extracted yield was mostly interpreted as on-resonance resonance 
thick target yield and translated to a resonance strength. There are large uncertainties 
associated with this approach, in particular with the determination of the resonance 
energies, which have not been unequivocally determined in the experiment. In some 
cases several of the quoted resonances agree with previously identified states, in 
other cases it needs to be confirmed that the observed yields really correspond to 
additional resonances and do not originate from tail contributions of resonant yield 
curves associated with the different states. As far as the resonance levels which have 
been observed in both studies are concerned the published strengths are compara- 
ble to each other. Nevertheless the reaction rate suggested by Vockenhuber et al. 
( Vockenhuber et al. 2007 ) is larger by more than a factor of two than the rates pro- 
jected on the basis of the in-beam gamma spectroscopy measurements, but it is in 
agreement with the projections by Rauscher et al. (Rauscher et al. 2000b[ l. The dif- 
ference is mainly due to the difference in resonance numbers. While the resonance 
identification in previous work was based on a careful analysis of the particular 
gamma decay characteristics of the observed levels, the analysis of the recoil data is 
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insufficient in providing information to differentiate between different resonances. 
It cannot be excluded that the number of identified states are overestimated; a more 
detailed gamma spectroscopy study with thin targets is therefore highly advisable to 
remove the existing uncertainties. 



9.5.2 Experiments with Neutron Beams 

A particularly interesting case is the origin of the long-lived gamma emitter 60 Fe. 
Its characteristic /-radioactivity has been observed with the INTEGRAL gamma ray 
telescope in supernova remnants near the solar system. These observations are com- 
plemented by recent AMS studies which su ggest high 60 Fe abundance in deep sea 
ferromanganese sediments ( jKnie et aT] |2004[ ). These 60 Fe observations have been 
interpreted as indication for the existence of a recent (« 3 million years) supernova 
event in the solar system vicinity. A more quantitative interpretation of the time 
and distance of the proposed supernova event requires a detailed knowledge of the 
nucleosynthesis history of 60 Fe. 

The radioactive 60 Fe isotope is produced by a sequence of neutron capture reac- 
tions of stable iron isotopes such as 58 Fe(n,y) 59 Fe(n,y) 60 Fe, the production rate and 
final abundance of the long-lived 60 Fe depends on the reaction rate of these feeding 
processes as well as on the rate of the 60 Fe(n,y) 61 depletion reaction. No experimen- 
tal information are available about the associated cross sections except for the neu- 
tron capture reaction 58 Fe(n,y) 59 Fe. Present simulations of the 60 Fe nucleosynthesis 
rely entirely of statistical model predictions of the neutron capture rates. Because 
of the relatively low level density in the associated 60 Fe, 61 Fe compound nuclei 
these model predictions are unreliable and need to be tested experimentally. This 
is underlined by the direct comparison between the experimental cross sections for 
neutron capture on the stable isotopes 56 Fe, 57 Fe, and 58 Fe which were all measured 
through neutron activation techniques and theoretical Hauser Feshbach predictions 
which show considerable discrepancies in particular in the cases of 56 Fe(n,y) 57 Fe 
and 57 Fe(n,y) 58 Fe. For 58 Fe(n,y) 59 Fe on the other hand , the agreement seems rea- 
sonable well but that cannot be extrapolated towards neutron captures on the more 
neutron rich Fe isotopes which are subject of the here proposed measurements. 

Particularly important is the determination of the reaction rate of 59 Fe(n,y) 60 Fe 
since it competes directly with the 59 Fe j3 -decay which would by-pass the produc- 
tion of 60 Fe. A direct measurement of this critical reaction in the traditional acti- 
vation or time of flight spectroscopy technique is not feasible because the target 
is radioactive and only small amounts can be accumulated, these small amounts 
nevertheless produce a large background activity level, which would prohibit any 
of the described methods. The cross section for the ground state decay branch of 
59 Fe(n,yo) 60 Fe can however be investigated using inverse 60 Fe(y,n) 59 Fe Coulomb 
dissociation techniques. The 60 Fe beam can be produced by fragmentation of a 
heavy ion such as 64 Ni on a light Be target at an energy of 500 MeV/u. The 59 Fe 
recoil products, and the released reaction neutrons, as well as yrays can be detected 
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with reasonable 200 keV resolution using a combination of a magnetic separator 
system and a neutron detector wall. This allows particle identification of all reaction 
products. 



9.5.3 Experiments with Radioactive Beams or Targets 



The depletion processes of long-lived radioactive isotopes includes the natural de- 
cay. Simulating this branch requires not only a good knowledge of the laboratory 
lifetime but also of the nature of the decay process since extreme environmental 
effects can change the decay rates drastically. In terms of j3 decay, the decay can 
be accelerated through the decay of thermally excited states as in the case of 26 Al. 
For decay through electron capture, the decay can be slowed down since the nuclei 
are completely ionized and the electrons have to be captured from the stellar plasma 
rather than from the inner K- or L-shell of the atom. This affects in particular the 
lifetime of 44 Ti, which primarily decays by electron capture. 

Often the depletion is primarily driven by nuclear reactions, such as 22 Na(p,y) 23 Mg, 
26 Al(p,y) 27 Si, 44 Ti(a,p) 47 V, or 60 Fe(n,y) 61 Fe, but also capture reactions on shorter- 
lived excited configurations of these nuclei are possible, such as 26 Al*(p,y) 27 Si. 
There are two possibilities for experimental studies of the reaction cross sections. 
The first one is based on the production of highly enriched long-lived radioactive 
targets, which can be prepared through standard chemical target preparation tech- 
niques using externally bred radioactive material, or by implantation of radioactive 
ions at low energy ISOL facilities. The disadvantage of both techniques is that the 
actual y measurements have to be performed in a high radiation background envi- 
ronment produced by the sample itself. 

Nevertheless, earlier measurements of reactio ns such as 22 Na(p,y) 23 Mg ( Seuthe et al. 



1990 1 and 26 Al(p,y) 27 Si ( Buchmann et al. 1984 1 relied entirely on this approach. In 



both cases a large number of resonances were detected and the resonance strength 
determined for calculating the reaction rates. The results for 22 Na(p,y) 23 Mg were 
confirmed by new direct measurements using improved target and detection tech- 
niques (Stegm uller et al.||1996[ ), resulting in the observation of an additional low 
energy resonance at lower energies. Complementary spectroscopy techni ques such 
as the study of the j3-delayed proton decay of 23 Al (Perajarvi et al. 2000) a nd the 



heavy ion reaction induced y decay of proton unbound states in 23 Mg (|Jenkins et al. 



2004 ) provided additional nuclear structure information which led to the reduction 
of uncertainties in the reaction rate. 

The situation is similar with 26 Al(p,y) 27 Si; after the initial study with radioac- 
tive targets (Buchmann et al. 1984[). A number of transfer experiments (Schmal- 



about the threshold levels in 21 Si not accessible to direct study by radiative cap- 
ture measurements lead to an improved reaction rate for 26 Al ground state capture. 
A first direct study of a lower energy resonance was successfully performed in in- 
verse kinematics at the ISAC facility at TRIUMF using the DRAGON recoil sepa- 



brock et al. 1986 Vogelaar et al. 1996| ) providing complementary information 
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rator ( |Ruiz et aT 2006)1 . The re sonance value is substantially smaller than the value 
quoted before (Vogelaar et al. 1996 1, which reduced the reaction rate slightly at 
temperatures anticipated for nova burning conditions. 

Not included in the reaction rate calculations are possible contribution s of proton 
capture on the thermally first excited state in 26 Al (Runkle et al. 2001 1. Recently 



number of indirect measurements have been performed to explore the possible con- 
tribution to the total reaction rate of 26 Al(p,y) 27 Si. Transfer reactions have been 
used to populate proton unbound states in 27 Si measuring the subsequent prot on de- 



cay to the ground state and the first excited state in 26 Al (Deibel et al. 2009). This 



approach allows to determine the branching and the relative strength of the proton 
decays for each of the unbound states. This can be used to scale the reaction rate 
component for the proton capture on the first excited state. 

Possible lo wer energy resonan ce contributions to the proton capture rates on the 
ground state ( |Lotay et aL 2009 1 and the excited state of 26 Al ( Lotay et al. 2009 1 
have been explored by y spectroscopy techniques probing the proton unbound state 
in 27 Si through heavy ion fusion evaporation reactions and measuring the y decay 
of proton unbound states. This is a particular efficient method to explore the levels 
near the threshold where proton decay is suppressed by the Coulomb barrier. The 
measurements provide critical information about spin and parity of the observed 
states but gives only limited information about the resonance strengths which is 
primarily determined by the proton decay strength. 

The main reaction for the depletion of 60 Fe in neutron rich environments is 
60 Fe(n,y) 61 Fe. The reaction rate used for nucleosynthesis simulations was for many 
years based on theoretical Hauser Feshbach model predictions. Recently an experi- 
ment has been performed at the FZK Karlsruhe in Germany to determine the stellar 
reaction cross section experimentally by neutron activation with the neutron beam 
resembling a quasi-stellar neutron spectrum (Uberseder et al. 2009). The activated 
60 Fe sample was prepared from PSI beamstop material. The cross section was de- 
termined from the characteristic 61 Fe y activity relative to the amount of 60 Fe nuclei 
in the target material. The latter was determined from the characteristic 60 Fe y ac- 
tivity of the target sample, based on this the experimental results suggest a cross 
section which is twice as large as standard Hauser Feshbach predictions suggesting 
a much more rapid depletion of 60 Fe in neutron rich environment than previously 
anticipated. The estimate of the number of 60 Fe nuclei, however relied on adopting 
a half-life of Ti / / 2 =l -49 Gy ( |Kutschera et al. 1984 ». Recent work suggested that the 
half-live is considerably larger model Ty^=lxil Gy ( |Rugel et"aL 2009 [ l. This would 
translate into a considerably larger amount of 60 Fe particles in the sample, suggest- 
ing a cross section which would be in fair agreement with the Hauser Feshbach 
predictions. New independent life time measurements for 60 Fe are clearly necessary 
to address this issue and remove the uncertainty in the interpretation of the radiative 
capture data. 
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